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I. INTRODUCTION 

Next year, the Large Hadron Collider (LHC) will begin 
operation at CERN, colliding protons at a center-of-mass 
energy of 14 TeV, seven times greater than that cur- 
rently available in pp collisions at Fermilab's Tevatron. 
The LHC luminosity should also be a factor of 10 to 100 
greater than the Tevatron's. The combined rise in energy 
and luminosity represents the opening of a new window 
into electroweak-scale physics. There will be copious pro- 
duction of heavy states such as electroweak vector bosons 
and top quarks, and, it is anticipated, Higgs boson(s) and 
physics beyond the Standard Model. 

Are we ready to exploit this new window? That is, 
is Standard Model physics at the LHC understood well 
enough to confidently extract physics beyond the Stan- 
dard Model? Of course the physics at the LHC is largely 
the physics of hard QCD processes (with important con- 
tributions from soft regimes as well). QCD itself will 
be probed in unprecedented regimes and with unprece- 
dented statistics. QCD governs the production of elec- 
troweak states in and beyond the Standard Model. Pure 
jet final states also need to be understood, because jet 
production rates are so large. Jets can fake electroweak 
signatures such as leptons or photons, and can fake miss- 
ing energy via mismeasurement. In general, QCD (plus 
electroweak) hard processes form significant backgrounds 
to almost all new physics searches at the LHC, as well as 
at the Tevatron. 

To be sure that our understanding of hard QCD at 
the LHC is good enough, it is important to check QCD 
predictions against current data from the Tevatron and 
HERA. The main experimental inputs are the strong cou- 
pling, a s (Mz), and the parton distribution functions. 

Recent progress in computing three-jet observables 
in e + e~ annihilation at next-to-next-to-leading order 
(NNLO) [3, to be described below, promises to improve 
the experimental uncertainty in a s as determined from 
e + e~ event shapes. On the other hand, the uncertainty 



in the current world average 



a 8 (Mz) =0.1182 ± 0.0027, 



(1.1) 



is already small, in comparison with other uncertainties, 
for almost all LHC processes. Parton distribution func- 
tions (pdfs), and their uncertainties, are critical to all 
QCD predictions, but as they were reviewed in the talks 
by Gwenlan Q and Diehl pj, I will not dwell on them 
here. Here I will focus more on our theoretical under- 
standing of the hard, short-distance structure of QCD 
processes, assessed when possible against Tevatron and 
HERA data for various processes and regimes. 

Because hard QCD is a vast subject, this talk will 
only scratch the surface. I will say next to nothing 
about several important subjects covered in part by 
other speakers, such as the high-energy (BFKL) limit, 
small- x physics and parton saturation [4| ; (hard) diffrac- 
tion 0, H[ ; and heavy quark production [3, @] • I will 
also not be able to cover multiple parton scattering and 
the underlying event; the substantial recent progress 
in matching leading-order QCD predictions with parton 
showers; other Monte Carlo developments; insights for 
collider physics being developed via soft collinear effec- 
tive theory; and new techniques for computing one-loop 
QCD amplitudes with many external legs. 

Section II of this article outlines the framework and ba- 
sic elements of fixed-order QCD computations and briefly 
discusses resummation. Section III describes a state- 
of-the-art application to e + e _ annihilation, the thrust 
distribution at NNLO [l|. In Section IV, developments 
in the theory of Higgs production at hadron colliders 
are summarized, along with selected decay channels and 
some backgrounds thereof. Section V is devoted to jet 
physics: definitions, substructure, rates and distribu- 
tions. Section VI covers the production of a vector boson 
in association with jets. Section VII discusses the pro- 
duction of a top quark pair plus a jet, and the tt forward- 
backward asymmetry. In Section VIII I conclude. 



II. PRELIMINARIES AND THEORETICAL 
TOOLS 

Asymptotic freedom in QCD Q guarantees that at 
short distances (large transverse momenta) the partons in 
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the proton are almost free, and are sampled essentially 
one at a time in hard collisions. This picture leads to 
the QCD-improved parton model, in which the hadronic 
cross section for production of a final state X factor- 
izes into products of pdfs f a and partonic cross sections 

-ab^X 



V" / dxidx 2 f a (xi,a s ,fj, F )f b (x2,a s ,n F ) 
» Jo 

(sxix 2 ;a s , hf, hr)- (2.1) 



Here \ip and fin are the factorization and renormaliza- 
tion scales, which are in principle arbitrary. In practice, 
truncating the cross section at a given order in perturba- 
tion theory induces dependence on fip and fj,R. 

Although parton distributions are nonperturbative 
quantities which must be measured experimentally at 
some short-distance scale /i, their evolution with /i is gov- 
erned by the DGLAP equation [8J, 



dfa{x,n) a s {p) f 1 df , 



whose kernel is known through NNLO 



(2.2) 



P, 



> ab {x,a s ) = P^+gpW^+g) P$(x)+0(a*). 

(2.3) 

The partonic cross section can be expanded similarly in 
powers of a s , 
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(a s ,fj,F,m) 



(0) , as(HR) ±.(i) r „ „ s 



2tt 



(2.4) 



where n a depends on the process. For typical collider 
processes, fin might be of order 100 GeV, for which 
a s {fiR) ~ 0.1. One might expect that the leading-order 
(LO), or Born level, terms in the expansion (ct(°)) would 
suffice to get a 10% uncertainty. However, for hadron col- 
lider cross sections, corrections from the next-to-leading 
order (NLO) terms in the a s expansion (a^) can in- 
crease the cross section by 30% to 80%. There are several 
reasons for the large corrections, some of which we shall 
discuss below. Thus, LO predictions are only qualitative; 
quantitative predictions require NLO corrections. If a 
few percent precision is desired, then the next-to-next- 
to-leading order (NNLO) terms, may also be required. 
Also one must be careful to describe the experimental 
setup (cuts, etc.) sufficiently accurately. 



A. Basic ingredients at fixed order 

What ingredients enter a perturbative QCD calcula- 
tion at LO, NLO, or NNLO? First of all, various partonic 
scattering amplitudes are required. These amplitudes are 
illustrated in fig. [T] for one of the simplest processes, the 
inclusive production of a Z boson at a hadron collider, 
followed by Z decay to an electron-positron pair. At LO, 
only tree amplitudes are needed. In this example, a sin- 
gle Feynman diagram contributes to qq — > Z — > e + e~. 
This diagram just needs to be squared, and convoluted 
with the pdfs, while incorporating any experimental cuts 
on the final state leptons. 

At NLO, one-loop amplitudes contribute to virtual cor- 
rections; for example, the one-loop correction to qq — > Z. 
The virtual corrections must be combined with real radi- 
ation; i.e., tree amplitudes having one additional parton 
in the final state. In the Z example, the subprocesses are 
qq — ► Zg, qg — » Zq and qg — ► Zq. The virtual and real 
corrections are separately divergent in the infrared (IR), 
which includes both soft and collinear regions. Usually 
the IR divergences are regulated dimensionally, by letting 
the number of spacetime dimensions be D = 4— 2e (with 
e < 0), and expanding both virtual and real contributions 
in a singular Laurent expansion around e = 0. There are 
1/e 2 singularities that cancel between virtual corrections 
and real corrections. Some of the 1/e singularities also 
cancel this way; others, representing initial-state collinear 
singularities, are absorbed into a renormalization of the 
pdfs. Ultraviolet poles are removed by coupling renor- 
malization. The finite remainder is then convoluted with 
the pdfs, as at LO. 

At NNLO, there are three types of terms: two-loop 
virtual corrections to the lowest-order process; mixed vir- 
tual/real corrections from one- loop amplitudes with one 
additional parton; and tree amplitudes with two addi- 
tional partons, as shown in fig. [TJ The IR cancellations 
are increasingly intricate, beginning now at order 1/e 4 . 

As the number of final-state partons in a process grows, 
the complexity of the theoretical computation increases, 
at every order in a s , but the issue is particularly prob- 
lematic at NLO and NNLO. Consider, as an example, the 
processes pp — > n jets. At LO, fast numerical programs 
allow the computation of the n-jet final-state up to about 
8 jets [l(3], depending on the computing time available. 
At NLO, there are no complete results for more than 
three jets. In this case, the main limitation is the lack of 
knowledge of the one-loop amplitudes for more than five 
external partons (except for one case, gg — * gggg). At 
NNLO, even the basic two-jet final state cannot yet be 
computed at NNLO. Here the required amplitudes are all 
known, and the main obstacle has been the integration 
over the singular final-state phase space. More interest- 
ing final states may include electroweak particles such as 
W, Z, or Higgs bosons in addition to jets. In each case, 
the state-of the-art value of n is the same or smaller as 
in the n-jet case, if one counts each electroweak particle 
as replacing one jet. 
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FIG. 1: Sample Feynman diagrams contributing to Z boson 
production at a hadron collider, at LO, NLO, and NNLO. 
Only one diagram is shown for each contributing amplitude, 
and some amplitudes are omitted. 




(a) (b) 

FIG. 2: Schematic depiction of two different types of subtrac- 
tion methods. Arrows represent momenta of hard external 
partons. (a) The dipole subtraction method is built around 
collinear singularities associated with individual partons. (b) 
The antenna subtraction method is built around soft singu- 
larities associated with pairs of partons. 



is finite, so it can be performed numerically, with 
generic experimental cuts. 



B. Singular phase-space integration 

Most of the recent advances in computing NNLO cor- 
rections to collider processes have come from the develop- 
ment of methods for integrating over singular regions of 
phase space where one or two partons are "unresolved" , 
i.e., are either soft or are collinear with another, hard 
parton. There have also been recent advances at NLO, 
in the context of matching parton showers with fixed- 
order computations, which exploit an understanding of 
the singular phase-space structure. 

At NLO, only one parton can be unresolved. For ex- 
ample, consider the final-state gluon in the tree-level am- 
plitude for qq — > Zg in fig. [TJ There is a soft singularity 
when the gluon momentum vanishes, k g — > 0. This sin- 
gularity can be interpreted classically, as radiation from 
the "accelerating" quark color charges when the quark 
and antiquark annihilate. Thus soft radiation is associ- 
ated with a pair of external partons. There are also two 
collinear singularities, one in which the gluon is collinear 
with the quark, from the Feynman diagram shown in 
fig. [U and one in which the gluon is collinear with the 
antiquark, from another Feynman diagram (not shown). 

Various techniques have been developed for handling 
these singular integrals at NLO, and some of them arc 
also being applied at NNLO. There are four general cat- 
egories of methods: 

• analytic — Usually carried out on a case-by-case 
basis, for the simplest processes only, with (at best) 
very simple cuts. 

• slicing [ll| — Thin strips are excised from the sin- 
gular regions of phase space. An approximate ver- 
sion of the cross section can be integrated analyti- 
cally in the strips. Outside the strips, the integral 



subtraction 0, [H, Q [H El E3 — One sub- 
tracts, over the entire phase space, a function mim- 
icking the exact cross section in singular regions, 
but which can be integrated more easily, even in 
the presence of cuts. The subtracted difference is 
integrated numerically. 



• direct numerical integration [18|, [19|, |2fJ — This 
approach can be carried out for a variety of pro- 
cesses, even at NNLO, but most effectively after 
using sector decomposition [2l[ to help separate the 
singularities. 

Within the category of subtraction methods, one can 
distinguish two different ways of building up subtraction 
terms for general processes out of simpler building blocks. 
The problem is that soft singularities connect all possible 
pairs of legs, while collinear singularities are associated 
with individual legs. Also, soft and collinear regions over- 
lap, so any smooth subtraction term should have a soft 
part and a collinear part. Figure HJa) sketches the build- 



ing block for the dipole subtraction method [151 ]. This 
function captures the collinear behavior near a particu- 
lar hard parton, plus part of the soft behavior connect- 
ing that parton to other color-correlated partons. 1 Fig- 
ure [2ljb) illustrates the building block for the antenna 
subtraction method [16j| . ft captures the soft behavior 
associated with a pair of color-connected partons, plus 
part of the collinear behavior near each of the members 
of the pair. 



1 The term "dipole" used here differs from the "dipole shower" 
used in the Monte Carlo community, which is closer to an antenna 
pattern. 
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C. Parton showers and NLO matching 

Parton showers represent an approximation to the 
soft and collinear radiation pattern, a resummation of 
leading logarithms (at present), which can be imple- 
mented probabilistically. Parton showers are a key part 
of Monte Carlo simulation programs such as PYTHIA [22j 
and HERWIG [23], which produce hadron-level events and 
are essential to experimental analyses. To improve the 
accuracy of parton-shower Monte Carlo programs, it is 
important to match them to fixed-order results. LO 
matching is an important subject which I cannot do jus- 
tice to here; see however, ref. [24j for a recent comparative 
study of different approaches. 

For even better accuracy, though, NLO matching is 
necessary. The program MC@NL0 [25j was the first to ac- 
complish NLO matching for a variety of different pro- 
cesses. It operates within the HERWIG Monte Carlo, which 
uses an angular-ordered shower. Implementing addi- 
tional processes in MC@NL0 has been nontrivial. An im- 
portant issue is to avoid the double-counting of emissions, 
between (i) the first step of the shower, and (ii) the ex- 
act radiation pattern in QCD, from which subtraction 
terms have been removed in the course of the NLO com- 
putation. If the shower radiation pattern differs from the 
form of the subtraction terms, a correction is necessary. 

Recently there have been some advances in achieving 
NLO matching in a more process-independent way. The 
POWHEG method [26], in contrast to HERWIG, generates the 
hardest parton emission first, and works well with p^- 
ordered showers. Recently it has been formulated for a 
general NLO subtraction method, and in par ticular for 
the Frixione-Kunszt-Signer [l4| and dipole [15[ methods, 
for generic processes [27j |. 

It has also been recognized [28j that subtraction meth- 
ods automatically supply radiation patterns that are cor- 
rect in the soft and collinear limits, and hence can be used 
to construct a parton shower. In this case the double- 
counting problem is solved automatically, because the 
first step of the shower coincides with the NLO subtrac- 
tion function. Two independent implementations of a 
shower based on dipole subtraction have appeared very 
recently [29|, as well as one (VINCIA) based on antenna 
subtraction [3fj]. With the recent flourishing of different 
parton-shower algorithms and NLO matching routines, 
it will be very interesting to compare their outputs for 
benchmark processes over the next year or two. 

D. Subtraction at NNLO 

Let us return now to the status of fixed-order re- 
sults. The NLO technique that has been carried out 
for the widest variety of collider processes is the dipole 
subtraction method. Programs such as MCFM [3l| and 
NL0JET++ |32j cover a variety of hadron collider, ep col- 
lider and e + e~ processes, limited mainly by the avail- 
ability of virtual corrections (one-loop amplitudes). The 



method has been generalized to handle massive final- 
state partons (such as top quarks) [33| : and a fully au- 
tomated version of the method in the massless case has 
appeared very recently [34j . 

It is natural to try to generalize this method to NNLO. 
Such a generalization is highly non-trivial, because there 
are now several different types of singular phase-space 
regions, for the contributions with two additional radi- 
ated partons, exemplified by the rightmost terms on the 
NNLO line in fig. [1] There can be singularities when 
both additional partons are soft, when one is soft and the 
other collinear; when both are collinear with a third par- 
ton; and when there are two independent collinear pairs. 
These different regions overlap with each other. Never- 
theless, progress has been made in constructing dipole- 
type subtraction terms for e + e~ annihilation to jets [35| . 
The method has also been adapted to give NNLO results 
for the inclusive production of a Higgs boson via gluon 
fusion at hadron colliders (36j |. 

The most complex process that has been treated to 
date at NNLO is that of three-jet event shapes in e + e~ 
annihilation, such as the thrust distribution [l[ to be de- 
scribed further in the next section. These results relied 
on building antenna-type subtractions to handle all the 
singular regions, and evaluating their phase-space inte- 
grals m Hi . 



E. Sector decomposition at NNLO 

Iterated sector decomposition is a strategy for doing 
singular integrals by partitioning the integration region, 
and then remapping it, in order to make the singularities 
one-dimensional. A very simple example of the procedure 
starts with the integral 

/= f f dxdy-^f^. (2.5) 
Jo Jo {x + yy 

Integrals like / are encountered in the NNLO corrections 
to vector boson production (for example) , from the inter- 
ference of two different initial-state radiation graphs in 
which a parton, radiated from either incoming line, splits 
into two additional partons. Although / is not singular 
as either x or y approaches zero with the other variable 
held fixed, there is a singularity as x = y — » 0. This 
singularity can be exposed by splitting the unit square 
into sectors A and B shown in fig. [3l and then remap- 
ping region A back to the unit square using x' = x, 
1 — y' = (l — y)/(l — x), and B back to the unit square us- 
ing y' = y, 1 — x' = (l—x)/(l — y). These transformations 
map the x — y — > singularity onto one variable, x' — ► 
or y' — ► 0, depending on the sector. This technique 
was first applied to multi-loop integrals [2l[ , and later to 
phase-space integrals [f8|, [ljl [2(| ■ Several iterations and 
many sectors may be required for state-of-the-art NNLO 
results. The expansion in e for one-dimensional singular- 
ities is straightforward, involving standard "plus" distri- 
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FIG. 3: Sector decomposition splits and remaps integration 
regions in order to expose multi-variable singularities. 



butions. Hence arbitrary observables can be integrated 
over phase space. 



F. Soft-gluon resummation 

In particular kinematic regions, fixed-order perturba- 
tion theory breaks down. This breakdown can be due to 
a mismatch between the kinematics of virtual and real 
corrections, enhanced by the strength of soft-gluon emis- 
sion. For example, in Z production by qq annihilation, 
the transverse momentum of the Z, q^(Z), vanishes for 
all virtual corrections, but is nonzero for real corrections, 
which are enhanced for small qT(Z). The leading behav- 
ior at L loops is 



da 



dq\{Z) 



(C F a s )- 



ln 2L-l „ 2 



3t(Z) 



q\(Z) 



(2.6) 



plus terms with fewer logarithms. Here Cf = 4/3 is the 
quark color charge. These large corrections imply that 
transverse-momentum resummation [381 ] is required for 
q T (Z) < M z . 

Distributions in the threshold variable z = M^/s q q 
behave similarly. Virtual corrections have z = 1, while 
real corrections have z < 1, and are enhanced as z — ► 1, 
with leading behavior, 



da 

dz 



Tn^-^l-z) 
1 - z 



(2.7) 



Such singular distributions can give rise to large correc- 
tions to inclusive production cross sections if they are 
convoluted with steeply falling parton distributions, ne- 
cessitating threshold resummation [39j . 

The kinematics are similar for production of the Higgs 
boson via gluon fusion, except that the color charge for 
gluons, Ca = 3, is much larger than for quarks, and 
the gluon pdfs are falling much faster than the quark 
pdfs in the relevant x range, x ~ 1CP 2 . Hence both 
and threshold resummation effects are much larger in this 
case than for vector boson production. 



one from e + e annihilation. The thrust 40] is a classic 
infrared-safe e + e~ event-shape observable, defined by 



T 



max Yh Wi 

h 



(3.1) 



The sum is over the final-state hadrons (or partons, in a 
perturbative calculation) with spatial momenta pi, and 
n is a unit vector, varied over all directions. There is a 
wealth of data to compare to, for yfs from 14 to 206 GeV, 
with the best statistics gathered on the Z pole. In the 
extraction of a s (Mz) from a fit of NLO QCD predictions 
to e + e~ event shapes, it is a long-standing problem (see 
e.g. ref. [4l[) that the central value depends on the ob- 
servable used, and on the range in that observable used 
in the fit. The error is dominated by the truncation of 
the perturbative series at NLO. 

Finally, 27 years after the first NLO event-shape re- 
sults Jl2| , the first NNLO results have appeared this 
year [Ij. At the Z pole, the effect of the NNLO correc- 
tions is to increase the NLO thrust distribution by about 
15-20% in the range 0.04 < 1 — T < 0.33. (The two-jet 
region, 1 — T < 0.03, requires resummation. The region 
1 — T > 0.33 cannot be produced by qqg final states, so 
it is also less perturbatively stable.) The relative uncer- 
tainty in the perturbative prediction is reduced by about 
30-40% with respect to NLO. 

The increase in the thrust distribution should lead to 
a somewhat smaller a s {Mz)\ but at NLO the thrust led 
to a larger than average value, a s (Mz) ~ 0.126. A new, 
more precise value of a s {Mz) requires experimental re- 
analysis, incorporating: 

1. a resummation of large logarithms for 1 — T — > 0, 



2. an analysis of power corrections of the 
Aqcd / Q using data off the Z pole, 



form 



3. other event-shape observables, which have also 
been computed (very recently) at NNLO [421 ]. 

Indeed, the first NNLO determination of a s from event- 
shape variables, using ALEPH data for six different vari- 
ables (including thrust) from several center-of-mass en- 
ergies, has just been reported (post LP07) [43| . 



a s (M z ) = 0.1240 ±0.0008 sta t± 0.0010, 

±0.0011 had ± 0.0029 thco . 



CXp 



(3.2) 



III. AN e+e" APPLICATION - THRUST AT 
NNLO 

Next I turn to some recent applications of these the- 
oretical techniques to collider processes, beginning with 



The result is a bit higher than the world average (jl.ip . 
The perturbative uncertainty, "theo" , has been cut 
roughly in half with respect to NLO. Further analyses, 
incorporating resummation and power corrections, are 
eagerly awaited! 
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IV. HIGGS PRODUCTION AT HADRON 
COLLIDERS 

A. Gluon fusion total cross section 

Let's begin the discussion of hadron collider applica- 
tions with Higgs boson production and decay. These re- 
sults are of great phenomenological importance for the 
Higgs search at the Tevatron, and particularly at the 
LHC. The gluon-fusion channel, gg — > H, via a top quark 
loop, represents one of the simplest final states, apart 
from vector boson production, so the theory has already 
been pushed to rather high order. At the same time, 
the process illustrates methods of analysis that should 
eventually be applied to more complex states. 

In the most likely mass range for the Standard Model 
Higgs boson, from 114 GeV to about 200 GeV, gluon 
fusion dominates the total production cross section an- 
It has been known since the early 1990s that the NLO 
corrections to <jh were huge, increasing it by roughly 
80% [4J] . Because the lowest-order process proceeds by a 
loop diagram, the NLO corrections already require two- 
loop integrals. To go to NNLO, the large-mt approxi- 
mation has been used. This approximation shrinks the 
top-quark triangle to a point, replacing it by an effective 
operator, HGf lv G> lva [45[, with coefficient C H - It re- 
duces the number of loops and the number of mass scales 
in the problem by one. This makes it feasible to perform 
all the NNLO computations, including the phase-space 
integrals, analytically for the case of inclusive produc- 
tion [46| . Threshold logarithms, of the form 



da 



H 



dz 



(c A a s y 



ln fc (l - z) 



1 



(4.1) 



with k = 0, 1, . . . 2L — 1, play a big role in the large pos- 
itive corrections. They have been resummed to next-to- 
next- to-leading logarithmic (NNLL) accuracy (47j. Nev- 
ertheless, a sizable uncertainty of order 10-15% remained 
in the prediction. 

Some progress has been made on reducing this uncer- 
tainty. As a spinoff from the computation of the NNLO 
DGLAP kernel for gluon evolution, Pjjf(x) 9], Moch 
and Vogt [H| were able to extract the leading singular 
terms (|4.1[) at three loops, L = 3 and k = 0, 1, 2, 3, 4, 5. 
They defined an approximation N 3 LO a p prox which is 
missing just the <5(1 — z) term, plus all terms that are 
nonsingular as z — > 1. They also employ the N 3 LO 
corrections to the coefficient Ch 0- The results are 
shown in fig. [5] The renormalization-scale dependence is 
now stabilized near the Higgs mass, with a residual un- 
certainty of order 5%. More recently, a similar analysis 



appr 



with the six 
. The results 



has been carried out at order "N 4 LO 
leading terms in eq. (|4.1[) at that order [5C 
change by about 2% with respect to the N a L0 a pp rO x re- 
sults shown in the figure. The total Higgs production 
cross section is now becoming a precision observable. 



o(pp -> 

-\ M H = 


, , 1 

H+X) [pb] 
120 GeV ^ 




NLO 


- N 3 LO appro 

N 2 LO 


LO ' : 



a(pp H> H+X) [pb] 
M u = 240 GeV 



0.2 0.5 
H r /M H 



2 3 



0.2 0.5 



2 3 



FIG. 4: Variation of the total Higgs production cross section 
with respect to the renormalization scale fi r , at various levels 
of approximation. Figure from ref. [4£|. 



B. Gluon fusion with decay to two photons 



Unfortunately, the total Higgs production cross sec- 
tion cannot be observed experimentally. For any final 
state, various experimental cuts have to be imposed. For 
example, for the two-photon final state from the decay 
H — > 77 there are cuts on the photons' transverse mo- 
mentum, rapidity, and isolation (with respect to hadronic 
energy). These cuts can now be mimicked at the parton 
level at NNLO, for gluon fusion production (in the large- 
mi limit) followed by H — > 77. (The same is now also 
true for the decay mode H -> WW -> Ivlv [M [H|.) 
Two independent programs are available for the 77 fi- 
nal state. The program FEHiP [20j] uses direct numerical 
integration, following sector decomposition for the real 
corrections. The more recent program HNNL0 (3(| em- 
ploys a subtraction method, based on the fact that at 
nonzero Higgs gx the problem is really a NLO calcula- 
tion; plus a knowledge, based on resummation studies, of 
the universal behavior of the gx distribution as q-j — ► 0. 

For regions of final-state phase space that are accessi- 
ble at LO, both programs show good perturbative sta- 
bility in going from NLO to NNLO. However, the behav- 
ior near boundaries of the LO-accessible region is more 
unstable, again reflecting the breakdown of fixed-order 
perturbation theory. One such boundary occurs at a 
large rapidity difference between the two photons [20I ]. 
At qt = 0, the photon pt cuts and fixed Higgs mass 
vtlh put a bound on how forward and backward the de- 
cay photons can go. More obvious boundaries occur in 
the transverse momenta of the two photons, pTmin and 
PTmax, which each have to be less than mff/2 at LO. One 
can see from fig. [5] that the pTmax distribution gets sig- 
nificantly stiffer, and the prmin distribution gets softer, 
from NLO to NNLO [II]. (The distributions are forced 
to be identical at LO.) 
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FIG. 5: Distributions of transverse momenta pt for the min- 
imum and maximum pt photons for Higgs production at the 
LHC, followed by the decay H -> 77, at LO, NLO and NNLO. 
Figure from ref. [36} |. 

C. The two-photon background 

In searching for a bump from the Higgs boson in the 
invariant-mass of photon pairs, it is also useful to know 
the characteristics of the continuum QCD background. 
The 77 background has many components, including: 

• electrons, positrons and hadrons faking photons 

• photons arising from copious ir° decays and elec- 
tron or positron bremsstrahlung 

• photons arising from fragmentation — radiation at 
small transverse momentum with respect to a jet 

• hard QCD radiation of photons. 

The first three categories of backgrounds can be sup- 
pressed fairly effectively by photon isolation cuts. 

Direct production of photon pairs in hard QCD begins 
at LO with the quark-annihilation process qq — > 77. The 
process qg — > 77(7 enters at NLO. It is heavily enhanced, 
by the large gluon pdf at small x, and by a final-state 
collinear singularity between each photon and the out- 
going quark. (At very small angles, the singularity is 
absorbed into the fragmentation contribution; it is also 
suppressed by isolation cuts.) At NNLO, the process 
99 ~ * 77: mediated by a virtual quark loop, enters for 
the first time. It is significant, however, because the 
gluon pdf enters twice. NLO programs for the inclu- 
sive di-photon background, including fragmentation ef- 
fects, were constructed in the late 1990's 52, 53]. The 
contribution from the gg — > 77 subprocess at its next-to- 
lcading order in the pdf-enhanced channel gg — ► 77(g) 
was added a few years later (54[. This year the quark- 
box contribution to qg — > 77(7 was included, as well as 
a resummation of the transverse momentum q-j of the 
di-photon pair at NNLL accuracy [H, [5rl. 

CDF has presented 207 pb _1 of data [57| on pairs of 
photons produced at the Tevatron with p T > 14 GeV 



pp > n*. *S - ' 96 TcV pp > ffX, \S - 1 .% lev 
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FIG. 6: Distributions in invariant mass Q and azimuthal angle 
separation A0, for pairs of photons produced at the Tevatron. 
Figure from ref. [561 ]. 



and di-photon invariant masses ranging from 10 to 100 
GeV. This range is useful for assessing the hard QCD 
background to the Higgs search at the LHC. If one scales 
a di-photon invariant mass of, say, 150 GeV, and typical 
photon pt cut of 40 GeV, down by the ratio of beam 
energies between the LHC and the Tevatron (seven), one 
gets a di-photon invariant mass of order 20 GeV, and pt 
cut of 6 GeV, not too far from the CDF cuts. So similar 
ranges of x values for the pdfs are being sampled. Fig- 
ure [6] shows the CDF data, in comparison with the NLO 
plus NNLL-resummed predictions [56[ . The agreement is 
generally quite good, within the available statistics, ex- 
cept for the region of small A<j>. Significant contributions 
to the small A<f) region can arise when a single parton 
yields both photons, e.g. gq — ► gq'y, with hard photon 
radiation off the final-state quark, followed by the frag- 
mentation of that quark to a second photon. These frag- 
mentation contributions are NLO (order a 2 a s ), and are 
not included in the computation of r ef. 1561] - They have 
been included in the program DIPH0X [53|, which does fit 
the CDF data at small A(j>. The small Acj> region will 
probably not be too important for the Higgs search; such 
large boosts are kinematically disfavored. 

In general, then, the di-photon background at the LHC 
seems to be in relatively good shape; although a com- 
putation of the qg — > 773 channel at its NLO, leading 
eventually to a computation of qq — > 77(33) at NNLO, 
would certainly be welcome. 



D. Vector boson fusion 

The second largest Higgs boson production mecha- 
nism, weak boson fusion (WBF), qq — + qqH, features 
a pair of forward tagging jets, as shown in fig. [3a). The 
NLO QCD corrections that dominate in the forward limit 
are sketched in fig. [TJb) . They only involve one quark line 
at a time, and have been known for a while to be quite 
modest, only 5% or so [58|,l59|. Recently, the complete set 
of QCD corrections to WBF were computed, as well as 
the one-loop electroweak corrections to all channels [60| . 



8 




H 





H 




(b) 




FIG. 7: (a) LO diagram for weak boson fusion, (b) Most 
important NLO QCD corrections, (c) Example of an addi- 
tional NLO QCD correction, (d) Sample NLO electroweak 
correction. 
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FIG. 8: Distribution in azimuthal separation A(f> for tagging 
jets in the gluon-fusion-plus-two-jets process, after applying 
WBF cuts, at LO (multiplied by the NLO K-factor), and at 
NLO. Figure from ref. [gj]. 



distribution is fairly stable, only flattening slightly. 



An example of an additional QCD correction is shown 
in fig. it can interfere with the graph in fig. [T^a) if 
the two quarks are identical. There are also s-channel 
annihilation graphs. However, these contributions are 
kinematically disfavored after typical WBF cuts (empha- 
sizing forward jets) are imposed. Ref. [601 ] confirms that 
these contributions are tiny. On the other hand, the one- 
loop electroweak corrections are found to be sizable and 
negative, of order —7%, so they must be included along 
with the QCD corrections. 

Higgs boson production via WBF has various back- 
grounds, of course, depending in general on the decay 
mode. It also has a source of "background" that is in- 
dependent of the decay mode, coming from the produc- 
tion of a Higgs boson via the gluon fusion subprocess 
(Hgg interaction), plus the radiation of at least two more 
jets. This "background" could in principle impact stud- 
ies of the WBF production mechanism. While most of 
the gluon-fusion-plus-two-jets background is eliminated 
by WBF cuts, it is important to understand how various 
distributions for this subprocess are affected by higher- 
order QCD corrections. Recently, gluon-fusion-plus-two- 
jets was computed at NLO, in the large-mt limit [6l| . 
The computation employed a semi-numerical evaluation 
of the one- loop virtual corrections [13] • The overall rate 
for this subprocess, with typical WBF cuts, increases by 
30% in going from LO to NLO. This increase is much 
less than that for the inclusive gluon-fusion process, but 
still significant. The azimuthal separation A<p of the two 
tagging jets is an incisive probe of the WBF production 
mechanism [59[ . Figure [8] shows how this distribution 
changes in going from LO to NLO, in the gluon-fusion- 
plus-two-jets "background" subprocess. The normalized 



V. JETS 
A. Jet definitions 

Jets are by far the most copious high-transverse- 
momentum objects produced at hadron colliders. They 
are common as well to almost all ep final states and a 
majority of e + e~ final states at high energy. Hence a 
thorough understanding of their production rates and 
properties in various regimes is highly desirable. 

There are two popular classes of algorithms for defin- 
ing jets. At e + e~ colliders, cluster algorithms have tra- 
ditionally been used [63| . A "distance" metric is defined 
between pairs of particles, which vanishes when the two 
particles are collinear or one is soft. The pair with the 
smallest distance is clustered into a proto-jet, and the 
process is iterated, until all proto-jets are separated by 
more than a specified distance, called the jet resolution 
parameter y. This "bottom up" procedure automatically 
assigns every particle to a jet. It is infrared-safe at the 
parton level: Adding an arbitrarily soft gluon, or split- 
ting one parton into two very collinear partons, does not 
change any of the subsequent clustering steps, so the final 
set of jets remains the same. This means that jet rates 
can be calculated (in principle) to any order in perturba- 
tion theory. 

Events at hadron colliders typically contain, in addi- 
tion to a hard partonic scattering process, an "underlying 
event" , consisting of forward beam remnants and parti- 
cles produced more centrally by interactions of specta- 
tor partons with each other and with partons from the 
hard process. The fcx algorithm 64] is a cluster-style al- 
gorithm that has been adapted for use at both ep and 
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FIG. 9: Fraction of hard events that fail an infrared-safety 
test, for various cone algorithms. Figure from ref. [68}. 



hadron colliders. It clusters particles that are "closer" to 
the beam axis into a beam jet, rather than into one of 
the other jets. However, some energy from the underly- 
ing event will still be associated with non-beam jets, and 
quantifying this is an issue, although it can be done [65} . 

Traditionally, hadron collider jet algorithms have been 
based instead on cones. Cone jets are circles of fixed ra- 
dius R in the plane of pseudorapidity rj and azimuthal 
angle <fr containing a set of final-state hadrons, or a set 
of calorimeter towers, with total transverse momentum 
p 3 ^ above a particular lower threshold. There are also 
rules for splitting and merging nearby cones in order to 
arrive at a stable configuration. Usually the cones are 
"seeded"; i.e., there is a prescription for specifying which 
circle centers will be considered first. The seeds might be 
towers with transverse momentum above a seed thresh- 
old. Minimum-bias events distribute energy uniformly in 
the r}-4> plane, making for a simple prescription for cor- 
recting for the underlying event in a cone algorithm. See 
e.g. refs. [H, [(^J for more details. 

Unfortunately, for seeded cones there is a danger of 
infrared unsafety: The final jet configuration can some- 
times change by a lot, with the addition of an arbitrarily 
soft gluon. Introducing a midpoint seed between all the 
previous seeds can reduce the problem, but not elimi- 
nate it. In the past year, a new, computationally practi- 
cal, seedless cone algorithm, SIScone, was invented [68j ] . 
This algorithm is infrared safe, as shown in fig. [5] The 
figure also indicates the fraction of hard events that fail 
an infrared-safety test, for some other popular cone algo- 
rithms. 



Packages are now available that give the user the flex- 
ibility to reconstruct jets with a variety of algorithms 
and parameters, and to do comparisons of jet properties 
between different algorithms [63, [H, [6^| . 

B. Jet substructure 

To what extent can the identities of underlying partons 
be deduced from properties of the jets they produce? 
In particular, can we distinguish light quark jets from 
gluon jets? It is not possible to do this event by event; 
however, it is possible to do it statistically, by study- 
ing how energy is distributed within the jet. In more 
detail, event kinematics are used to select gluon-rich or 
gluon-depleted samples of jets. For these samples, one 
can measure the distributions of the jet-shape function 
*i>(r/R) — the fraction of energy for a jet with cone size 
R that is found in a smaller cone with r < R. If the 
distributions differ for the gluon-rich and gluon-depleted 
samples, then one could use the jet-shape distributions 
to separate quark and gluon jets, statistically. (Similar 
studies have been carried out in the past using sub-jet 
multiplicities, in e + e~ [j(J, ep [7l| andpp [72j] collisions.) 

In pp collisions, CDF studied the dependence of 
the jet-shape distribution on the transverse momen- 
tum of jets reconstructed with the midpoint algorithm 
(R = 0.7) [73]. Higher transverse-momentum jets are 
gluon-depleted, in comparison with lower transverse- 
momentum jets. Good agreement was found with 
PYTHIA [22j, for transverse momenta ranging from 
50 GeV (73% gluon, 27% quark, according to PYTHIA; 
(*(0.3/i?)) = 0.705 ± 0.015) all the way up to 350 GeV 
(20% gluon, 80% quark; (*(0.3/i?)) = 0.93 ±0.02). 

Recently, ZEUS studied the jet-shape distribution in 
ep collisions, with 368 pb _1 of data, using two-jet events 
and taking the lowest transverse-momentum jet in order 
to select a gluon-enriched sample [74| . Figure [10] shows 
the ZEUS data, using the fcx algorithm with resolution 
parameter D = 2.5 and jet transverse momenta in the 
narrow range between 14 and 17 GeV. The agreement 
with NLO predictions (from the program DISENTjl5|) is 
excellent, to within a percent or so. Thus the substruc- 
ture of jets in generic events in both pp and ep collisions 
is understood very well now. 



C. Inclusive jet rates 

The production rate for very high transverse- 
momentum jets at the Tevatron provides a direct test of 
QCD interactions at the shortest possible distance scales. 
Of course, the rate also depends on the parton distri- 
bution functions, in particular the large- a: gluon distri- 
bution. Because jet production rates are rapidly-falling 
functions of pr, they are very sensitive to the jet en- 
ergy scale calibration, as well as to tails in the jet en- 
ergy resolution. Both D0 and CDF have presented new 
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FIG. 12: Ratio of inclusive-jet data from CDF to NLO theory 
for the midpoint cone algorithm (R — 0.7), plotted versus pt 
for various bins in rapidity Y. Figure from ref. [77| . 



FIG. 10: Jet-shape distribution in ep collisions for jets in one- 
jet events, vs. the lowest pr jet in two-jet events, at the same 
jet pt- Figure from ref. [741 ]. 
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measurements of inclusive-jet rates using midpoint cone 
algorithms with a cone radius of R = 0.7, and based on 
approximately 1 fb -1 of data. 

Figure [TT1 shows data from DO divided by theory, for 



central jets (|yj e t| < 0.4) and binned in px- The theory 
includes NLO, plus an estimate of the NNLO terms based 
on threshold logarithms c omp uted at next-to-leading log- 
arithmic (NLL) accuracy |76| . NLO theory without these 
corrections is about 10% lower in this pt range. The cen- 
tral value of the cross section from DO is closer to the 
threshold-enhanced theory; however, the systematic un- 
certainty is too large to distinguish the two curves. Fig- 
ure [T3] shows that the analogous data from CDF [73] is 
also in good agreement with NLO theory [32j], at both 
central and forward rapidities. 

For both the DO and CDF results, the systematic un- 
certainty is dominated by the jet energy scale uncertainty 
and exceeds the statistical error for essentially all points. 
However, the uncertainty arising from the pdfs is even 
larger, and increases with pt- This fact illustrates that 
the high-px jet data provide a strong constraint on pdfs, 
and can be used to reduce the pdf uncertainties in par- 
ticular regimes, particularly the gluon distribution g{x) 
at large x. 



D. Di-jet azimuthal distributions 

At leading order in QCD, two jets produced in a pp 
collision should emerge back-to-back in the plane trans- 
verse to the beam axis, i.e. with the maximum azimuthal 
angle between them, A(f> = it. QCD radiation pushes the 
azimuthal angle to smaller values. In a multi-jet event, 
the azimuthal angle is defined to be between the two 
highest pt jets. For A(f> « n, multiple soft gluon ra- 
diation dominates, and fixed-order perturbation theory 
cannot be trusted. Because of the way the azimuthal 
angle is defined, three parton final-states can only pro- 
duce Acf> > 2tt/3. The A(f> distribution is an excellent 
test of how well QCD describes complex hadronic final 
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FIG. 14: Azimuthal decorrelations from HI in the forward 
region. Figure from ref. p79( ] . 



VI. VECTOR BOSONS PLUS JETS 

The production of a single vector boson in associ- 
ation with multiple jets is a background to searches 
for supersymmetry at the LHC. For example, the final 
state Z + n jets, when the Z decays to vv, is a back- 
ground to the missing-transverse-momentum plus multi- 
jet searches. CDF and DO have studied such events at 
the Tevatron, for the decays W — > tv and Z — > £ + £~, 
plus up to 4 jets in the final state [8l|. Recently, CDF 
presented data [82] on Z plus 1, 2 and 3 jets, in com- 
parison with fixed-order LO and NLO predictions. Fig- 
ure [T5] shows the excellent absolute agreement, to within 
10%, with NLO predictions for 1 and 2 jets, the maxi- 
mum number for which NLO results are currently avail- 
able [83J. These results highlight the importance of ex- 
tending NLO theory to larger numbers of jets, in order to 
be ready for the enormous data sets that will be available 
in these channels at the LHC. 



states, because it is relatively insensitive to the overall 
jet energy scale and to the pdfs. 

Figure [TBI shows D0's measurement of the azimuthal 
decorrelation in di-jet events at the Tevatron [781 ] . com- 
pared to LO and NLO predictions from NL0JET++ [32] . 
LO theory fails for both small and large values of A<j>, 
for the reasons mentioned above. However, NLO theory 
does very well for Acf> < 2ir/3 (even though it is effec- 
tively an LO calculation in this regime), and also pushes 
the agreement for large values considerably closer to tt. 

In contrast, NLO theory is unable to describe the az- 
imuthal correlations measured in ep collisions by HI [7{|, 
shown in fig. [T4j at least for small values of Xsj- The 
NLO three-jet predictions [8(| are closer to the data than 
are the two-jet predictions, but they still do not produce 
enough decorrelation at small xbj- Does this signal a 
breakdown of fixed-order perturbation theory for these 
kinematics, and a need for small- x resummation? Per- 
haps. On the other hand, HI also compared the data 
with a few different models incorporating small- a; resum- 
mations, and all such models were found to be too low 
for small A0 as well. 



VII. TOP QUARKS PLUS JETS 

The final-state tt+jet is another important background 
to supersymmetry at the LHC. The cross section is large, 
and the additional jet can boost the ti system so that neu- 
trinos from top quark decay generate large missing trans- 
verse momentum. The NLO corrections to this process 
were computed recently [84| . They require the evaluation 
of many virtual Feynman diagrams, including pentagon 
diagrams with massive propagators, and a large number 
of subtraction terms for the real corrections [33j. The 
results greatly reduce the uncertainty on the tt+jet cross 
section at the LHC. 

At the Tevatron, a pp collider, the forward-backward 
asymmetry of ti pairs is an interesting observable that 
probes the dynamics of top quark production. It can be 
defined by 

A t = N t (y t > 0) - N t (y t < 0) 

Afb N t (y t >0) + N t (y t <0)> ^ 

where Nt(yt > 0) is the number of top quarks produced 
with positive rapidity. This quantity vanishes at leading 
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this conference. The observable used by CDF [6j and 
D0 [a, |88| is a bit different from eq. (|7.1|) . It employs 
the rapidity difference between the t and t, 



A t = N t(Vt >Vt)-N t {y t < yi ) 
N t (y t > Vt) + N t (y t < y { ) ' 



(7.2) 



The difference between these two observables for the 
NLO inclusive definition was studied recently [87[. For 
the same choice of pdfs, cuts, etc., it is found that 
A FB = (5.1 ±0.6)%, whereas A 1 = (7.8 ±0.9)%. The ob- 
servable A* is somewhat larger than A FB because events 
where both the t and t go forward can contribute to A* 
but not to ^4p B - 



CDF measures, with 1.7 fb of data, 
A 1 = (28±13 st at±5 syst )%; 
while D0 measures, with 0.9 fb -1 of data, 



4' 



(12 ± ; 



± 1 



syst)% , 



(7.3) 
(7.4) 



the latter number has not been corrected for reconstruc- 
tion effects. The statistical errors are still large, of course. 
However, in light of the large central values, it will cer- 
tainly be interesting to follow these results as more data 
are analyzed. 



FIG. 15: Differential cross sections for Z/j*(—+ e + e ) + 1 or 
2 jets at the Tevatron, as a function of p ^*, a nd in comparison 
with NLO theory JH- Figure from ref. [§g]. 

TABLE I: tt forward-backward asymmetry at the Tevatron, 
from refs. [H|l|. 



^fb(%) 


tt inclusive 


tt+jet inclusive 


ttOj exclusive 


LO 





-6.9 





NLO 


3.8 


-1.5 ±1.5 


6.4 



order, 0{a 2 s ), just like the forward-backward asymmetry 
for muons in the QED process e + e~ — > ji + a" . But it 
is nonvanishing at NLO, or 0{a 3 s ) [H, [H, , as shown 
in table HI It is also nonvanishing for the process tt±jet, 
computed at its leading order, again 0{cP s ). Here a jet is 
required to have p^* > 20 GeV. Ref. computed the 
asymmetry (|7.1[) in tf±jet events at NLO, O(a^), and 
found the striking result that the asymmetry is drasti- 
cally reduced, from —6.9% to essentially zero. One might 
wonder what this result portends for the NNLO value of 
the tt inclusive asymmetry — will it too receive a large 
correction? 

In the mean time, the first forward-backward asym- 
metry measurements from the Tevatron were reported at 



VIII. CONCLUSIONS 

This talk has surveyed some of the recent progress over 
the past year or so in our quantitative theoretical under- 
standing of hard QCD processes at colliders. An increas- 
ing number of processes are now known at NLO, and 
for a few benchmark processes, NNLO precision is avail- 
able. Experiments at the Tevatron and at HERA test 
such QCD predictions over a wide range of kinematics. 
The experiment between experiment and theory is gen- 
erally very good, except near kinematic boundaries such 
as small transverse momentum and small x, for which re- 
summations and reorganizations of the perturbation the- 
ory should be performed. It is critical to push the "loops 
and legs" frontier to bring more processes to NLO and 
NNLO accuracy, and also to incorporate such processes 
into parton showers, while retaining that accuracy. But 
in general, we are "virtually" ready for the startup of the 
LHC next year! 
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